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Abstract: An analysis is presented of events containing jets including at least one b- 
tagged jet, sizeable missing transverse momentum, and at least two leptons including a pair 
of the same electric charge, with the scalar sum of the jet and lepton transverse momenta 
being large. A data sample with an integrated luminosity of 20.3 fb“^ of pp collisions 
at y/s = 8 TeV recorded by the ATLAS detector at the Large Hadron Collider is used. 
Standard Model processes rarely produce these final states, but there are several models 
of physics beyond the Standard Model that predict an enhanced rate of production of 
such events; the ones considered here are production of vector-like quarks, enhanced four- 
top-quark production, pair production of chiral 6'-quarks, and production of two positively 
charged top quarks. Eleven signal regions are defined; subsets of these regions are combined 
when searching for each class of models. In the three signal regions primarily sensitive 
to positively charged top quark pair production, the data yield is consistent with the 
background expectation. There are more data events than expected from background in 
the set of eight signal regions defined for searching for vector-like quarks and chiral b'- 
quarks, but the significance of the discrepancy is less than two standard deviations. The 
discrepancy reaches 2.5 standard deviations in the set of five signal regions defined for 
searching for four-top-quark production. The results are used to set 95% CL limits on 
various models. 
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1 Introduction 

The Standard Model (SM) has been repeatedly confirmed experimentally. Nonetheless 
there is a need for physics beyond the SM (BSM) at about the TeV scale, with additional 
features that explain the baryon asymmetry of the universe, specify the nature of dark 
matter, and provide a mechanism to naturally stabilize the Higgs boson mass at its observed 
value of approximately 125 GeV [1, 2]. This paper reports on a search for BSM physics 
resulting in pairs of isolated high-transverse-momentum (high-px) leptons^ with the same 
electric charge, hereafter denoted as same-sign leptons, (or three or more leptons of any 
charge) missing transverse momentum, and 6-jets. This is a promising search channel since 
the SM yields of such events are small, and several types of BSM physics may contribute. 

Among the models that predict enhanced same-sign lepton production are those that 
postulate the existence of vector-like quarks, an enhancement of the four-top-quark pro¬ 
duction cross section, the existence of a fourth generation of chiral quarks, or production 
of two positively charged top quarks. A common data sample is used to search for each of 
these signatures, but separate final selection criteria are defined based on the characteristics 
of each signal model. 

^Only electrons and muons are considered in the search. Tan leptons are not explicitly reconstructed, 
but electrons and muons from r decay may enter the selected samples. 


- 1 - 






B 

T 

Mass (model) 

Wt 

Zb 

Hb 

Wb 

Zt 

Ht 

0.50 TeV (singlet) 

42 

31 

27 

50 

17 

33 

0.50 TeV (doublet) 

100 

0 

0 

0 

34 

66 

0.55 TeV (singlet) 

43 

30 

27 

49 

18 

32 

0.55 TeV (doublet) 

100 

0 

0 

0 

37 

63 

0.60 TeV (singlet) 

44 

29 

26 

49 

19 

31 

0.60 TeV (doublet) 

100 

0 

0 

0 

38 

62 

0.65 TeV (singlet) 

45 

29 

26 

49 

20 

30 

0.65 TeV (doublet) 

100 

0 

0 

0 

40 

60 


Table 1: B and T quark branching fractions (in percent), assuming the singlet and (T, B) 
doublet models of ref. [25]. In the doublet case it is assumed that the mixing of the T 
quark with the Standard Model bottom quark is much smaller than the mixing of the B 
quark with the top quark. 


Several extensions to the SM that regulate the Higgs boson mass in a natural way 
require the existence of vector-like quarks (VLQ) [3-21], where ‘vector-like’ means that the 
left- and right-handed components transform identically under the SU(2 )l weak isospin 
gauge symmetry. Since quarks with this structure do not require a Yukawa coupling to the 
Higgs field to attain mass, their existence would not enhance the Higgs boson production 
cross section, and thus the motivation persists for a direct search [22]. There are several 
possible varieties of VLQ; those having the same electric charge as the SM b- and t-quarks 
are called B and T. In addition the exotic charge states T 5/3 and B_^i^ may occur, where 
the subscripts indicate the electric charge. Vector-like quarks may exist as isospin singlets, 
doublets, or triplets. Arguments based on naturalness suggest that VLQ may not interact 
strongly with light SM quarks [23, 24]. Thus it is assumed for this analysis that VLQ 
decay predominantly to third-generation SM quarks. For the B and T quarks, charged- 
and neutral-current decays may both occur {B —>■ Wt,Zb, or Hb] T —)• Wb,Zt, or Ht), 
providing many paths for same-sign lepton production for events with BB or TT pairs. 

The branching fractions to each allowed final state are model-dependent, and the ones 
occurring in models where the B and T exist as singlets or as a (T, B) doublet [25] are used 
as a reference. These branching fractions vary with the B or T mass, and values for some 
masses are given in table 1. Since the pair production of heavy quarks is mediated by the 
strong interaction, the cross section is identical for vector-like quarks and b' quarks (de¬ 
scribed below) of a given mass. The next-to-next-to-leading-order (NNLO) cross sections 
from TOP-1 — h v2.0 [26, 27] are used in this paper. The T 5/3 quark must decay to W^t, 
and therefore both single and pair production of this quark can result in same-sign lepton 
pairs, and both sources are considered. 

Same-sign lepton pairs may also arise from the production of four top quarks (titt). 
The SM rate for this production is small (~ 1 fb [28, 29]), but there are several BSM physics 


- 2 - 









models that can enhance the rate, such as top compositeness models [30-32] or Randall- 
Sundrum models with SM fields in the bulk [33]. These can generically be described in 
terms of a four-fermion contact interaction with coupling strength C' 4 t/A^, where Ci^t is 
the coupling constant and A is the scale of the BSM physics [32]. The Lagrangian for this 
interaction is 


(^rT^^r) (^R7/i^R) (1-1) 

where is the right handed top spinor and the 7 ^ are the Dirac matrices. Two specific 
models are also considered. The first is sgluon pair production, where sgluons are colour- 
adjoint scalars that appear in several extensions to the SM [34-39]. If the sgluon mass is 
above the top quark pair-production threshold, the dominant decay is to tt, resulting in four 
top quarks in the final state^ {tttt). The cross sections considered in this paper are rescaled 
to the next-to-leading order (NLO) prediction of Ref. [41]. The second model is one with 
two universal extra dimensions under the real projective plane geometry (2UED/RPP) [42]. 
The compactification of the extra dimensions leads to discretization of the momenta along 
their directions. The model is parameterized by the radii R 4 and R 5 of the extra dimensions 
or, equivalently, by tukk = 1/^4 and ^ = R 4 /R 5 . This model predicts the pair production 
of tier^ (1)1) Kaluza-Klein excitations of the photon with mass approximately 

V^rriKK that decay to ti with a branching fraction assumed to be 100%. The model also 
predicts a four-top-quark signal from tiers (2,0) and ( 0 , 2 ). Cosmological observations 
constrain mxK in this model to lie approximately between 600 GeV and 1200 GeV [43]. 

A fourth generation of SM-like quarks includes a charge —1/3 quark, called the b' [44- 
47]. Under the assumption that the 6 '-quark decays predominantly to Wt, b' pair produc¬ 
tion results in four W bosons in the final state. If two W bosons with the same electric 
charge decay leptonically, there will be a same-sign lepton pair in the final state. If the 
6 '-quark can also decay to Wq, where g is a light (u or c) quark, some b' pairs would also 
result in same-sign lepton pairs or trileptons (provided that at least one ft'-q^ark decays 
to Wt), and therefore the possibility of such decays is explored as well. The existence of 
additional chiral quark generations greatly enhances the Higgs boson production cross sec¬ 
tion, so if the new boson observed at the LHG is a manifestation of a minimal Higgs sector, 
additional quark generations are ruled out [48-53]. However, a more complex Higgs sector, 
as in some Two-Higgs-Doublet models [54], allows a fourth generation of chiral quarks. 

Production of two positively charged top quarks via uu —)• tt can also result in an excess 
of same-sign lepton pairs. This process may be mediated via s- or t-channel exchange of a 
heavy particle [55, 56]. In the t-channel exchange case, the process must include a vertex 
with a flavour-changing neutral current (FGNG). The neutral particle that is exchanged 
may be a vector, Z-like, particle or a scalar, Higgs-like, particle. Past searches for a new 
Z' boson have already put strong constraints on this possibility, thus only the scalar case 
is considered, with the following generic model Lagrangian [57]: 

^FCNC = HutHtHu + HctniHc -\- h.c. (1.2) 

^The decays predominantly to ti is model-dependent, as discussed in ref. [40]. 

® A tier of the Kaluza-Klein towers is labelled by two integers, corresponding to the two extra dimensions. 
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Figure 1: Leading-order diagrams for (a) vector-like top quark pair production, (b) sgluon 
pair production, and (c) same-sign top quark pair production through a BSM flavour¬ 
changing Higgs coupling. 


where H is a Higgs-like particle with mass niH and HutH and HctH denote the flavour¬ 
changing couplings of H to up-type quarks. Two scenarios are tested, one corresponding 
to a possible FCNC coupling of the newly discovered Higgs boson {mu = 125 GeV) and 
the other to a second scalar boson with a mass in the range [250, 750] GeV. If the mass 
of the mediating particle is much greater than the electroweak symmetry breaking scale, 
an effective four-fermion contact interaction can describe the process, thus extending the 
search to non-scalar particles. The corresponding Lagrangian contains separate operators 
for the different initial-state chiralities [58]: 

At = 

- ^^(■Ul7^A)Ar7Ar) “ ^^(«La7''Ab)AR67ARa) + h.C. 

where Gll, and Crr are the coefficients of effective operators corresponding to each 
chirality conhguration and A is the scale of the BSM physics. The Clr and C'-^^ terms 
lead to kinematically equivalent events, hence only one term is considered in this paper. 

Leading-order Feynman diagrams for the production in pp collisions of some of the 
signals searched for in this analysis are presented in hgure 1. 

Previous searches by the ATLAS collaboration [58] using an integrated luminosity of 
1.04 fb“^ of pp collisions at a centre-of-mass energy ^/s = 7 TeV and the GMS collabora¬ 
tion [59], using an integrated luminosity of 19.5 fb“^ of pp collisions at yfs = 8 TeV, did 
not observe a significant excess of same-sign dilepton production. The ATLAS result was 
used to set a lower limit of 450 GeV on the mass of a heavy down-type quark, under the 
assumption that the branching ratio to Wt is 100%, while the CMS result set upper limits 
on the four-top-quark production cross section of 49 fb,^ on the sum of the tt and tt pro¬ 
duction cross sections of 720 fb, and on the tt production cross section of 370 fb. The CMS 
collaboration used the same-sign lepton signature to search for quarks [60], ruling out 
such quarks with mass below 0.80 TeV, and as part of a broader search for vector-like T 
quarks [61], ruling out such quarks with mass less than 0.69 TeV. A more recent search by 

^Unless stated otherwise, all limits reported are at the 95% confidence level (CL). 
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the ATLAS collaboration [62] using an integrated luminosity of 20.3 fb“^ of pp collisions at 
y/s = 8 TeV with similar final states to those reported here was interpreted in the context 
of supersymmetric models. The present analysis improves upon the y/s = 7 TeV ATLAS 
analysis by using a larger data set recorded at a higher centre-of-mass energy, having a 
higher signal acceptance, and expanding the range of BSM models considered. 

2 Data and Monte Carlo simulations 

The data were recorded by the ATLAS detector [63] in LHC pp collisions at y/s = 8 TeV 
between April and December 2012, corresponding to an integrated luminosity of 20.3 fb“^. 
The ATLAS detector consists of an inner tracking system surrounded by a superconducting 
solenoid that provides a 2 T magnetic field, electromagnetic (EM) and hadronic calorimet¬ 
ers, and a muon spectrometer. The inner detector provides tracking information from pixel 
and silicon microstrip detectors within pseudorapidity^ \r]\ < 2.5, and from a transition ra¬ 
diation tracker that covers \r]\ < 2.0. The EM sampling calorimeter uses lead as absorber 
and liquid argon (LAr) as the active medium, and is divided into a barrel region that covers 
[t^I < 1.475 and endcap regions that cover 1.375 < |r^| < 3.2. The hadronic calorimeter 
consists of either LAr or scintillator tile as the active medium, and either steel, copper, or 
tungsten as the absorber, and covers \r]\ < 4.9. The muon spectrometer covers \rj\ < 2.7, 
and uses multiple layers of high-precision tracking chambers to measure the deflection of 
muons as they traverse a toroidal field of approximately 0.5 (1.0) T in the central (endcap) 
regions of the detector. A three-level trigger system selects events to be recorded for offline 
analysis. 

Signal and the background sources that contain prompt same-sign leptons or trileptons 
are modelled using Monte Carlo (MC) simulations. The remaining background sources are 
determined from the data, as described in section 4. B and T pair production is modelled 
using the protos v2.2 [25] generator using the MSTW2008LO [64] parton distribution 
functions (PDEs), with pythia v6.4 [65] used to model extra gluon emission and hadroniz- 
ation. T 5/3 production (both single and pair) is modelled with madgraph v5.1 [ 66 ] using 
the CTEQ6L1 [67] PDEs, with pythia v8.1 [ 68 ] used for hadronization. Production of 
four top quarks is modelled under four scenarios: i) Standard Model, ii) contact interac¬ 
tion, iii) sgluon pair, and iv) 2UED/RPP. The sgluon case is generated with pythia v6.4 
using the CTEQ6L1 PDEs; the other three models are generated with madgraph using 
the MSTW2008LO PDEs followed by pythia v8.1; in the case of 2UED/RPP the bridge 
generator [69] is used to decay the pair-produced excitations from madgraph to ti. The 
simulated 2UED/RPP samples correspond to the tier (1,1) for the symmetric {R 4 = R 5 ) 
case, with mxK ranging from 600 to 1200 GeV. Constraints on the asymmetric {R 4 > R 5 ) 

®ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) 
in the centre of the detector and the 2 -axis coinciding with the axis of the beam pipe. The x-axis points 
from the IP to the centre of the LHC ring, and the y-axis points upward. Cylindrical coordinates {r,(j>) are 
used in the transverse plane, (j> being the azimuthal angle around the beam pipe. The pseudorapidity is 
defined in terms of the polar angle 0 as ?7 = — lntan(0/2). For the purpose of the fiducial selection, this 
is calculated relative to the geometric centre of the detector; otherwise, it is relative to the reconstructed 
primary vertex of each event. 
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case are derived by an extrapolation that uses kinematical considerations [70]. These con¬ 
siderations also permit the extrapolation to signals arising from tiers (2,0) and (0,2) from 
the generated tier (1,1) signal. Pair production of b' events is modelled with the pythia 
v8.1 generator for b' masses ranging from 400 to 1000 GeV, using the MSTW2008LO PDFs. 
Production of two positively charged top quarks via a contact interaction is also modelled 
using PROTOS [71] and pythia v6.4, with three different chirality configurations of the 
contact interaction operator; production via the FCNC exchange of a Higgs-like particle is 
modelled with madgraph with pythia v8.1 used for showering and hadronization. The 
MSTW2008LO PDFs are used in simulating both types of tt production. 

The background contributions from ttW and ttZ (abbreviated as ttW/Z hereafter) 
and ttW~^W~ are modelled with madgraph followed by pythia v6.4, while WZ and 
ZZ plus jet production and W^W^jj production are modelled using sherpa vl.4 [72]. 
Background from the production of three vector bosons is modelled using madgraph and 
PYTHIA v6.4, and backgrounds from ttH, WH and ZH production are modelled using 
PYTHIA v8.1. The CTEQ6L1 PDFs are used for the ttWjZ, three-vector-boson, WH and 
ZH samples, the CTIO [73] PDFs are used for the WZ, ZZ, W^W^jj and ttH samples, 
and the MSTW2008LO PDFs are used for the tiW^W~ sample. In most cases (excluding 
background contributions that are negligibly small) the cross sections are scaled to match 
next-to-leading-order calculations. 

A variable number of additional pp interactions are overlaid on simulated events to 
model the effect of multiple collisions during a single bunch crossing, and also the effect of 
the detector response to collisions from bunch crossings before or after the one containing 
the hard interaction. Events are then weighted to reproduce the distribution of the number 
of collisions per bunch crossing observed in data. The detector response is modelled using 
either a geant4 [74, 75] simulation of the entire detector or a geant4 simulation of the 
inner tracker and of the muon spectrometer combined with a fast simulation of shower 
development in the calorimeter [76]. Some samples are generated with both types of sim¬ 
ulation, to allow direct comparison between the two, and agreement was found within the 
systematic uncertainty assigned to the efficiency estimate. In all cases the simulated events 
were reconstructed using the same algorithms that were applied to the collision data. 

3 Event selection 

The final states considered in this search require the presence of two leptons with the same 
electric charge in the event (events with additional leptons beyond the same-sign pair are 
also accepted). In addition, two or more jets are required, at least one of which is consistent 
with origination from a 6-quark, and sizeable missing transverse momentum E™'®® is also 
required, indicating the presence of neutrinos coming from W boson decays. The criteria 
used for each of these objects are given below. 

Each event is required to pass either an electron trigger (where the chosen triggers 
require either an isolated electron with pT > 24 GeV or an electron with pT > 60 GeV with 
no isolation requirement) or a muon trigger (where the triggers chosen require either an 
isolated muon with px > 24 GeV or a muon with px > 36 GeV with no isolation require- 
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ment). The trigger efficiency for electrons is ~ 95% while for muons it is ~ 75%, resulting 
in trigger efficiencies that range from ~ 95% for events with two muons to > 99% for events 
with two electrons. In addition, events are required to have at least one reconstructed ver¬ 
tex, which must be formed from at least hve tracks with px > 0.4 GeV. If multiple vertices 
are reconstructed, the vertex with the largest sum of the squared transverse momenta of 
its associated tracks is taken as the primary vertex. Since the events used in this analysis 
tend to have vertices with many associated tracks, the correct vertex is selected in more 
than 99% of the events. 

Electrons are identified by requiring a track to match an electromagnetic calorimeter 
energy cluster, subject to several criteria on the shape of the shower and the consistency 
between the shower and track. The selection requirements are varied with the r] and pT 
of the electron candidate to optimize the signal efficiency and background rejection [77]. 
The track is required to be within 2 mm in z of the reconstructed primary vertex of the 
event. A hit in the innermost layer of the inner detector is required to reject photon 
conversions. Energy clusters in the calorimeter associated with an electron are required 
to have transverse energy > 25 GeV and jr/j < 2.47, with the barrel/endcap transition 
region 1.37 < |??| < 1.52 excluded. The candidate is required to be isolated from additional 
tracks within a cone of variable AR = (Ar/)^ -|- = 10 GeV/px [78], such that the 

sum of the transverse momenta of the tracks within that cone must be less than 5% of the 
electron px- In addition, electrons are required to be separated from any jet by at least 
AR = 0.4. 

Muons [79] are identified from hits in the muon system matched to a central track, 
where the track must be within 2 mm in z of the primary vertex, and are required to have 
an impact parameter in the transverse plane that differs from the beam position by less 
than three impact parameter standard deviations. Requirements are placed on the number 
of hits in various layers of the muon system, and on the maximum number of layers where 
hits are missing. Muon pairs that are consistent with the passage of a cosmic ray are 
discarded. Muons are subject to the same track-based isolation requirement as electrons. 
Muons are also required to be separated from any jet by AR = 0.04 -|- 10 GeV/px, and to 
have Px > 25 GeV and \rj\ < 2.5. Events with a muon within A(f> x A6 = 0.005 x 0.005 of 
any electron are rejected. At least one of the selected leptons is required to match a lepton 
identified by the trigger. 

Jets are reconstructed from energy clusters in the calorimeter using an anti-fej al¬ 
gorithm [80-82] with radius parameter 0.4. If one or more jets are within AR = 0.2 of 
an electron, the jet closest to the electron is discarded (i.e. the cluster of energy in the 
calorimeter is treated as an electron rather than a jet). To suppress jets that do not ori¬ 
ginate from the primary vertex in the event, the jet vertex fraction (JVF) is defined by 
considering all tracks with px > 0.5 GeV within the jet, and finding the fraction of the 
summed px from tracks that originate from the primary vertex. Jets with px < 50 GeV 
and \ri\ < 2.4 that are matched to at least one track are required to have JVE greater than 
0.5. All jets are required to have px greater than 25 GeV (after energy calibration [83]) 
and jpj less than 2.5. 

A multivariate algorithm [84] is used to test if each jet is consistent with having arisen 
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from a 6-quark, based on the properties of the tracks associated with the jet. A requirement 
is placed on the output of the discriminant such that « 70% of 6-quark jets and « 1% of 
light-quark or gluon jets pass in inclusive simulated ti events. All jets that meet this 
criterion are called ‘6-tagged’ jets. 

The missing transverse momentum is calculated as the negative of the vector 
sum from all calorimeter energy clusters, with jet and electron energy calibrations applied 
to clusters associated with those objects, and corrected for the energy carried away by 
identified muons. Energy scale corrections applied to electrons and jets are also propagated 
to Events are required to have E™'®® > 40 GeV. 

If the same-sign leptons are both electrons, their invariant mass rriee is required to be 
greater than 15 GeV and to satisfy \mee — rnz{= 91 GeV)| > 10 GeV. These requirements 
reject events from known resonances decaying to an electron-positron pair where the charge 
of either the electron or positron is misidentified. Finally, the scalar sum of all jet and 
lepton transverse momenta is required to be greater than 400 GeV, since the signals 

considered here produce a high number of particles with high transverse momenta. These 
preselection criteria are applied to all searches; some of them are tightened when optimizing 
the selection for each signal model (see section 5). Figure 2 shows the distributions of 
and E™'®® after applying this selection (except for the requirements on and E™'®® 
themselves). 
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Figure 2: Distributions of the missing transverse momentum and scalar sum of jet 

and lepton transverse momenta ffT after applying the preselection criteria, for events with 
one, two, or more than two 6-tagged jets. The points with error bars are the data, the 
stacked histograms show the expectation from background, and the shaded band is the 
total uncertainty on the background expectation. 
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4 Background estimation 


Background arises from two distinct sources: SM processes that result in same-sign lepton 
pairs, and instrumental backgrounds where objects are misidentified or misreconstructed 
such that events appear to have the required set of leptons. The former category includes 
production of W^W^jj, tiW/Z, tiW+W-, tiH, WH, ZH, tWZ, tH, kkZ and ZZ with 
a heavy-flavour jet, or three vector bosons. In addition, the four-top-quark production 
predicted in the SM is included as a background to all searches for other signals, though 
its contribution is small due to the small cross section. All of these processes have small 
cross sections, and their expected yields are computed using simulation. Known differences 
between the lepton selection and 6-tag efficiencies between the MC simulation and data 
control samples are taken into account when computing the expected yields. 

Instrumental backgrounds have contributions from two categories; i) events where 
one or more jets are misidentified as leptons, or which contain non-prompt leptons, and ii) 
events that contain two leptons of opposite charge, where one of the charges is mismeasured. 
The ‘matrix method’ is used to estimate the contribution from events with misidentified 
(fake or non-prompt) leptons. In this method, the default (‘tight’) lepton identification 
criteria (section 3) are relaxed to form a ‘loose’ sample. Lepton isolation requirements are 
not imposed, and therefore the loose sample contains a larger fraction of fake/non-prompt 
leptons than the tight sample. The fraction of real leptons (meaning prompt leptons from 
the decay of a W, Z, or H boson) passing the loose criteria that also pass the tight criteria 
is referred to as r. Similarly, the fraction of fake/non-prompt leptons passing the loose 
cuts that also pass the tight cuts is referred to as /. Using measured values of r and /, one 
can construct a matrix that relates the observed yields of dilepton events in the categories 
loose-loose, loose-tight, and tight-tight to the real-real, real-fake, and fake-fake yields. 
An analogous procedure is applied to three lepton events starting with categories for all 
possible combinations of three loose or tight leptons, resulting in an estimate of the number 
of events with one or more misidentified leptons in the selected sample (Aj^^j^g). 

Single-lepton events are used to measure r and /. The criteria used to select these 
events are different for electrons and muons due to the differences in the sources of fake/non¬ 
prompt leptons for each flavour. For electrons, r is measured using events with K™'®® 
> 150 GeV, where the dominant contribution is from W —>■ eu, and / is measured using 
events with the transverse mass of the and electron® mT(IU) < 20 GeV and -|- 

mT(VF) < 60 GeV, where the dominant contribution is from multijet production where 
one or more jets is misidentified as an electron. For muons, r is measured using events with 
mT(VF) > 100 GeV, a sample dominated hyW—^ fiu, and / is measured using events where 
the impact parameter of the muon with respect to the primary vertex is more than five 
standard deviations from zero, consistent with muons arising from heavy-flavour hadron 
decays. The small contribution of real leptons to the control samples used to measure / 
is estimated from simulation, and this contribution is subtracted from the sample. The 

®The transverse mass of a lepton and is defined as rriTiW) = — cos A</)) where pTt. is 

the lepton px and A</) is the azimuthal angle between the lepton and the direction of the missing transverse 
momentum vector. 
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values of r and / are parameterized with respect to properties of the leptons (e.g. |??| and 
Pt) and of the event (e.g. the number of 6-tagged jets). Typical values are r = 0.90 and 
/ = 0.20 - 0.40 for electrons, and r = 0.95- 1.00 and / = 0.12 -0.30 for muons. 

The triggers used for low-pT leptons require isolation; since the tight and loose lepton 
criteria differ in their isolation requirements, fake/non-prompt leptons in events where only 
the low-pT triggers fired are more isolated, on average, than those from an unbiased trigger, 
meaning that / for these leptons is substantially higher. Therefore, r and / are measured 
separately for samples collected with the different triggers, and the appropriate values are 
applied based on the lepton triggers that fired in each event. A further complication may 
arise due to the small number of events in the loose sample, which can lead to a calculated 
value of that is negative or very close to zero. In the case of negative values 

is set to zero when computing limits. To properly estimate the statistical uncertainty 
on the fake/non-prompt lepton contribution given the small number of events, a Poisson 
likelihood for the estimate from the matrix method is used, and the standard deviation of 
the probability density function (p.d.f.) from this likelihood is used to set the uncertainty. 
In cases where the prediction from the matrix method is less than or near zero, the standard 
deviation is computed relative to zero rather than to the mean of the p.d.f. 

Charge misidentification (‘Q mis-Id’) is negligible for muons due to the small probab¬ 
ility for muons to radiate photons, the long lever arm to the muon system and the fact 
that the charge is measured in both the inner detector and the muon spectrometer. For 
electrons, the rate of charge misidentification is calculated from a sample of Z ^ ee events, 
selected with no requirement placed on the charge of the two electron tracks. It is assumed 
that the rate at which the charge of an electron is misidentified varies with the \r]\ and 
Pt of each electron but is uncorrelated between the two electrons in each event. Further 
assuming that the sample consists entirely of opposite-sign electron pairs, the number of 
measured same-sign events Nsi where one electron is in the {\v\^Pt) bin and the other 
in the bin is expected to be 


N^^N^^{e, + ej) (4.1) 

where iV®-^ is the total number of events in the i-j \r]\-pT bin, and e is the rate of charge 
mismeasurement. The value of e in each |??|-pT bin is then extracted by maximizing the 
Poisson likelihood for the observed number of same-sign pairs in each |7/|-pT bin to be 
consistent with the expectation from equation 4.1. One limitation of this estimate is that 
electrons from Z decay only rarely have large pT, rendering the uncertainty on the charge 
misidentification rate for high-px electrons large. To reduce this uncertainty, the rate of 
charge misidentification is estimated using simulated tt events as a function of the electron 
Pt- This rate is scaled to match the rate observed in data for the px range covered by 
the Z events, and the rate for electrons with larger px is extrapolated according to the 
scaled prediction from simulation. Closure tests comparing the number of events in the 
same-sign Z peak to the expectation based on the opposite-sign Z peak and the charge 
mismeasurement rates were performed in data and simulation and show good agreement. 
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To determine the number of events expected from charge mismeasurement in the sig¬ 
nal region, a sample is selected using the same criteria as for the analysis selection, except 
that an opposite-sign rather than same-sign ee or e/x pair is required. The measured e 
values are then applied to each electron in this sample to determine the expected number 
of mismeasured same-sign events in the analysis sample. One source of charge mismeasure¬ 
ment is from ‘trident’ electrons, where the electron emits a hard photon that subsequently 
produces an electron-positron pair, resulting in three tracks with small spatial separation. 
If the wrong track is matched to the EM cluster, the charge may be incorrect. However, 
such electrons would also appear to be isolated far less frequently than electrons that do 
not emit hard radiation. Therefore the value of r for trident electrons is lower than for 
electrons that have a correctly measured charge, meaning that they also contribute to the 
fake/non-prompt electron estimate from the matrix method. To avoid double-counting 
events with trident electrons in the background estimate, the charge mismeasurement rate 
is measured in a data sample where the non-prompt/fake contribution, estimated using the 
matrix method, has been removed. 

Simulation was used to estimate the sources of events in the signal regions that have 
fake/non-prompt leptons and/or electrons with mismeasured charge, and it is found that 
ti events provide the dominant contribution. 

The background estimates are validated using samples where one or more of the 
preselection criteria are vetoed so that the samples are statistically independent and the ex¬ 
pected yield from signal events is small. One such validation region, called the ‘low Hic+W 
region, is defined by applying the preselection criteria, except that the requirement on Hic 
is modified to 100 GeV < N't < 400 GeV. This validation region is particularly useful be¬ 
cause the background composition is similar to that of the preselection (including the fact 
that tt events are the dominant source of the fake/non-prompt lepton and charge mismeas¬ 
urement background contributions). The predicted and observed yields in this validation 
region are given in tables 2 and 3. Events with three leptons are considered explicitly in 
table 3 since the fake/non-prompt lepton background contribution from trilepton events 
is not negligible, and it is important to check that this component of the background is 
well understood. The “other bkg.” category includes WWW, WWZ, WH, ZH, ttWW, 
SM four-top-quark, and single-top-quark production. Similar agreement between the data 
yield and background expectation is observed in validation regions where no requirement 
on 6-tagged jets is imposed, and where is required to be less than 40 GeV or Et is 

required to be in the range 100-400 GeV. 
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Sample 

ee 

efj, 


Q mis-Id 

136 ± 2 ± 41 

118± 1±35 

— 

Fake/Non-prompt 

153 ± 11 ± 107 

225 ±11 ±158 

29 ± 3 ± 20 

tiWjZ 

4.57 ± 0.19 ± 1.88 

14.2 ±0.3 ±5.8 

8.43 ±0.27 ±3.56 

tlH 

0.39 ±0.04 ±0.04 

1.31 ±0.08 ±0.13 

0.76 ±0.06 ±0.07 

Dibosons 

5.57 ±0.45 ±1.08 

15.9 ±0.8 ±2.9 

9.00 ± 0.58 ± 1.79 

Other bkg. 

0.32 ±0.11 ±0.11 

0.75 ±0.20 ±0.20 

0.27 ±0.06 ±0.06 

Total bkg. 

299 ± 11 ± 115 

375 ± 11 ± 162 

47 ± 3 ± 20 

Data 

271 

307 

52 


Table 2: Observed and expected numbers of events in the validation region 

for the same-sign dilepton channels. The first uncertainty is statistical and the second is 
systematic (systematic uncertainties are described in section 6). 


Sample 

eee 

ee/x 

Fake/Non-prompt 

8.0 ±2.3 ±5.6 

13.2 ±2.4 ±9.2 

tiW/Z 

1.20 ±0.09 ±0.46 

2.55 ±0.13 ±0.87 

tlH 

0.07 ±0.02 ±0.01 

0.28 ±0.03 ±0.03 

Dibosons 

5.78 ±0.51 ± 1.14 

6.78 ±0.57 ±1.33 

Other bkg. 

0.04 ±0.02 ±0.02 

0.11 ±0.02 ±0.02 

Total bkg. 

15.1 ±2.4 ±5.7 

22.9 ±2.5 ±9.4 

Data 

15 

18 


Sample 

efifi 


Fake/Non-prompt 

17.9 ± 2.8 ± 12.5 

1.34 ±0.55 ±0.94 

tiWiZ 

3.38 ± 0.16 ± 1.15 

2.70 ±0.14 ±1.00 

tlH 

0.32 ±0.03 ±0.03 

0.14 ±0.02 ±0.01 

Dibosons 

8.42 ± 0.57 ± 1.78 

9.23 ±0.65 ±1.82 

Other bkg. 

0.12 ±0.02 ±0.02 

0.15 ±0.03 ±0.03 

Total bkg. 

30.1 ±2.8 ±12.7 

13.6 ±0.9 ±2.4 

Data 

36 

14 


Table 3: Observed and expected numbers of events in the low-HT^+lb validation region 
for the trilepton channels. The first uncertainty is statistical and the second is systematic 
(systematic uncertainties are described in section 6). 

5 Selection optimization 

The selection is defined to optimize the expected limit on signals. Since many BSM physics 
models (each of them dependent on mass and/or coupling parameters) could result in 
anomalous production of the sort sought in this analysis, defining a selection that is sensitive 
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to all of them is a challenge. As a first step toward a solution, the same-sign top signal is 
considered separately from the others, as it has unique characteristics: contributions are 
expected dominantly from positively charged lepton pairs^, and the jet multiplicity tends 
to be lower. The selection for same-sign top events is optimized with respect to ffx; ill™®®, 
and the number of 6-tagged jets N^,, with contributions from the ee, e^, and fifi channels 
considered separately. 

The remaining signals share a similar final-state topology, but the distribution of events 
differs between them in several variables. Therefore several event categories are defined, 
based on features of the events such as i/x, £1™®®, and N^,, as shown in table 4. Splitting 
the sample in this manner provides good overall efficiency for signal events, while allowing 
events that are least likely to arise from background (i.e. events with large values of i^x, 
or Nh) to be treated separately in the analysis, thereby enhancing the sensitivity to 
BSM physics. The boundaries between categories in ffx and were chosen to optimize 
the sensitivity to four-top-quark signals; these values are close to optimal for the other 
signals considered as well. 

All of the categories are considered when searching for vector-like quarks or chiral b'- 
quarks, while only the categories that require at least two 6-tagged jets are considered when 
searching for the production of four top quarks. One consequence of defining several signal 
categories is that the data-driven background estimates are subject to large statistical 
fluctuations. To mitigate this, all lepton flavours are summed within each category. The 
signal regions are defined based on the expected yields of signal and background, taking 
into account statistical and systematic effects, without considering the distribution of data. 


^Production of uu —>■ tt is also possible through these processes but the production cross section in pp 
collisions for this process is two orders of magnitude lower than that for uu —>■ tt. Therefore, considering 
only the tt final state reduces the background by a factor of two while having only a small impact on the 
signal. 
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Definition 


Name 


+ eee + ee^ + e/i^ + Nj > 2 



Nb = l 


SRVLQO 


400 <Ht < 700 GeV 

Nb = 2 

> 40 GeV 

SRVLQl 

SR4tO 


Nb> 3 


SRVLQ2 

SR4tl 


Nb = l 

40 < < 100 GeV 

SRVLQ3 



> 100 GeV 

SRVLQ4 


Ht > 700 GeV 

Nb = 2 

40 < Sip*"" < 100 GeV 

SRVLQ5 

SR4t2 


S“*"" > 100 GeV 

SRVLQ6 

SR4t3 


Nb >3 

> 40 GeV 

SRVLQ7 

SR4t4 


e+e+, e+//+, /i+//+, Nj G [2,4], A(/>f£ > 2.5 


> 450 GeV 


Nb > 1 


TT’miSS 

£/r-p 


>40 GeV 


SRttee, SRtte/i, SRtt/x/x 


Table 4: Definitions of the different signal regions. Nj is the number of jets that pass the 
selection requirements, and is the separation in cf) between the leptons. In regions 

SRVLQ0-SRVLQ7, contributions from all lepton flavours are summed. 

6 Systematic uncertainties 

Tables 5 and 6 show the sources of systematic uncertainties that contribute more than 1% 
uncertainty on the expected background or signal yield for the four-top/6YVLQ selection. 
These uncertainties have similar impact on the expected yields for the other signal models. 
For the yields derived from simulation, the largest source of uncertainty is the cross-section 
calculation. For the ttWjZ background, this is based on variations in the PDFs, variations 
of the renormalization and factorization mass scales (varied up and down by a factor of four 
from the nominal value of 172.5 GeV) [85], and variations in the parameters controlling 
the initial-state radiation model, resulting in a 43% uncertainty. For other background 
contrbutions, varying the renormalization and factorization scales results in uncertainties 
of 30% for WZ and ZZ production, 25% for W^W^jj production, -|-38%/ — 26% for 
tiW+W- production, and 10% for ttH, tH, WH, ZH, tWZ, WWW and ZIVVF produc¬ 
tion. These uncertainties, applied to the event yields shown in tables 8 and 9, result in the 
overall cross section uncertainties reported in table 5. The uncertainty on the integrated 
luminosity is 2.8% [86]. This uncertainty applies only to the backgrounds estimated from 
simulation, not to the data-driven estimates of the fake/non-prompt lepton and electron 
charge mismeasurement backgrounds, so the overall contribution of the luminosity uncer¬ 
tainty shown in table 5 is less than 2.8%. The largest detector-specific uncertainties arise 
from the jet energy scale [83], the 6-tagging efficiency [84], and the lepton identification 
efficiency [77, 79]. 

Systematic uncertainties on the background contributions estimated from data are 
evaluated separately. Six effects are considered when assigning the systematic uncertainty 
on the predicted yield of events from electron charge mismeasurement: i) the statistical 
uncertainty on the probability for an electron to have its charge mismeasured, ii) the 
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Source 

VLQ signal region number 

0 

1 

2 

3 

4 

5 

6 

7 

Cross section 

±8.0 

±13.6 

±15.1 

±11.1 

±12.1 

±16.8 

±25.2 

±23.8 

Jet energy scale 

+ 1.7 
-1.6 

+1.2 

-1.8 

+ 1.4: 
-1.7 

+1.8 

-2.1 

-1-2.6 

-4.2 

-H3.8 

-1.5 

-H8.5 

-4.8 

-H7.3 

-2.9 

6 -tagging efficiency 

±1.0 

±2.6 

+5.7 

-5.5 

-HI.9 
-2.0 

-HI.6 
-1.7 

-H3.8 

-3.7 

-H5.1 

-5.0 

-H8.3 

-8.2 

Lepton ID efficiency 

±1.3 

±1.6 

±1.6 

-H2.1 

-2.0 

-H2.1 

-2.0 

-H2.2 

-2.1 

-H2.8 

-2.2 

±2.5 

Jet energy resolution 

±0.5 

±0.2 

±3.1 

±1.9 

±0.3 

±0.9 

±0.8 

±3.4 

Luminosity 

±0.9 

±1.1 

±1.3 

±1.4 

±1.5 

±1.5 

±2.1 

±1.9 

Fake/non-prompt leptons 

±33 

±18 

±25 

±23 

±26 

±16 

±1.5 

±3.8 

Charge misID 

-1-5.9 

-5.7 

-1-9.3 

-9.1 

-1-5.4 

-5.1 

-1-7.4 

-6.7 

-H5.0 

-4.6 

-H8.7 

-8.1 

-H9.0 

-8.5 

-Hll.O 

-10.1 


Table 5: The largest systematic uncertainties (in %) on the total background yield for the 
four-top/6yVLQ selection. 


Source 

VLQ signal region number 

0 

1 

2 

3 

4 

5 

6 

7 

Jet energy scale 

-H11.3 

-9.0 

+11.5 

-6.3 

+28.0 

-17.3 

+3.7 

-2.1 

+5.4 

-2.4 

+3.9 

-2.0 

+4.5 

-6.5 

+6.6 

-3.0 

6 -tagging efficiency 

-H2.5 

-3.0 

+6.3 

-6.1 

+16.4 

-15.9 

+3.1 

-3.7 

+3.4 

-4.0 

+7.4 

-7.2 

+7.6 

-7.4 

+12.1 

-11.9 

Lepton ID efficiency 

±2.9 

±2.9 

±2.8 

±2.9 

+3.2 

-3.1 

±2.9 

+3.2 

-3.1 

+3.0 

-2.9 

Jet energy resolution 

±0.8 

±2.5 

±3.9 

±0.3 

±0.7 

±0.7 

±1.0 

±0.1 

Luminosity 

±2.8 

±2.8 

±2.8 

±2.8 

±2.8 

±2.8 

±2.8 

±2.8 


Table 6: The largest systematic uncertainties (in %) on the yield of a representative signal 
(600 GeV vector-like B pair production) for the four-top/6YVLQ selection. 


statistical uncertainty on the px-dependent scale factor, iii) the difference observed in sim¬ 
ulated Z boson events between the true charge mismeasurement rate and the rate obtained 
by applying the same method as is used for the data, iv) the difference in the px-dependent 
scale factor when measured using different ti simulated samples, v) the variation in the 
result observed when the width of the Z peak region is varied, and vi) the statistical uncer¬ 
tainty on the correction for the overlap in the measurement of charge misidentification and 
fake-electron background estimates. The magnitudes of these effects depend on the event 
characteristics, so the uncertainty on the background from electron charge misidentification 
varies from 23 to 40% in the signal and control regions, as presented in Tables 2 and 7-9. 
The expected yield of fake/non-prompt leptons is subject to uncertainties in the real and 
fake/non-prompt lepton efficiencies that arise from i) variations in the values of r and / 
when different control regions are used to measure them, ii) the small number of events in 
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those control regions, and in) the MC model used to subtract the real lepton contribution 
from the fake/non-prompt lepton control region. When assessing effect i, the following 
alternative control regions are used: for electrons, the alternative fake/non-prompt control 
region requires one loose electron and < 20 GeV, while for muons, the alternative 

control region requires one loose muon, m'^iW) < 20 GeV and ill™®® + mT(lV) < 60 GeV. 
In both cases the expected contribution from real leptons in the control region is subtrac¬ 
ted using simulation. The alternative control regions for r are formed by increasing the 
requirement on S™'®® from > 150 GeV to > 175 GeV for electrons and by increasing the 
requirement on mT(lV) from > 100 GeV to > 110 GeV for muons. Effects i)-iii) sum to 
a 70% uncertainty on the predicted yield of fake/non-prompt leptons. 

7 Results 

The observed yields for each signal selection are given in tables 7-9 and figure 3. The CLg 
method [87, 88 ] is used to assess the consistency between the observed yields and each 
potential BSM physics signal, where the log-likelihood ratio Lr is used as the test statistic. 
For each model, Lr is defined as 

LR = -21ogl^ (7.1) 

where Lg+b (Lb) is the Poisson likelihood to observe the data under the signal-plus- 
background (background-only) hypothesis. Pseudo-experiments are generated under each 
hypothesis, taking into account statistical fluctuations of the total predictions according 
to Poisson statistics, as well as Gaussian fluctuations in the signal and background ex¬ 
pectations describing the effect of systematic uncertainties. The quantities CLg^b and 



SRttee 

SRttep 

SRtt/ip 

ttw/z 

0.58 ± 0.06 ± 0.25 

1.20 ± 0.09 ± 0.53 

0.64 ± 0.07 ± 0.28 

tlH 

0.05 ± 0.02 ± 0.01 

0.12 ± 0.02 ± 0.02 

0.03 ± 0.01 ± 0.01 

Dibosons 

0.27 ± 0.14 ± 0.07 

0.38 ± 0.09 ± 0.10 

0.19 ± 0.12 ± 0.04 

Fake/Non-prompt 
Q mis-Id 

Other bkg. 

0.87 ± 0.79 ± 0.61 

2.66 ± 0.25 

0.01 ± 0.08 ± 0.00 

2.92 ± 1.27 ± 2.04 

2.79 ± 0.26 

0.05 ± 0.08 ± 0.01 

0.34 ± 0.29 ± 0.24 

0.12 ± 0.11 ± 0.03 

Total bkg. 

4.5 ± 0.8 

7.5 ± 1.3 ± 2.5 

1.3 ± 0.3 ± 0.4 

Data 

6 

5 

2 

p-value 

0.38 

0.84 

0.45 


Table 7: Observed and expected numbers of events with statistical (first) and systematic 
(second) uncertainties for the positively charged top pair signal selection. The p-values for 
agreement between the observed yield and the expected background in each signal region 
are reported. 
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SRVLQO 


SRVLQl/SR4tO 

SRVLQ2/SR4tl 

ttwjz 

16.2 

± 

0.3 

± 

7.0 

12.6 

± 

0.3 

± 

5.4 

1.24 

± 

0.09 

± 

0.53 

tlH 

2.5 

± 

0.1 

± 

0.3 

1.8 

± 

0.1 

± 

0.2 

0.26 

± 

0.03 

± 

0.05 

Dibosons 

11.2 

± 

0.6 

± 

2.8 

0.95 

± 

0.19 

± 

0.25 

0.07 

± 

0.12 

± 

0.05 

Fake/Non-prompt 

42.1 

± 

5.4 

± 

24.6 

8.61 

± 

2.34 

± 

5.02 

1.17 

± 

0.82 

± 

0.68 

Q mis-Id 

20.8 

± 

0.7 

± 

5.2 

15.1 

± 

0.6 

± 

3.5 

0.74 

± 

0.11 

± 

0.18 

Other bkg. 

1.76 

± 

0.13 

± 

0.17 

0.75 

± 

0.04 

± 

0.10 

0.10 

± 

0.08 

± 

0.03 

Total bkg. 

94.5 

± 

5.4 

± 

24.9 

40.0 

± 

2.4 

± 

7.3 

3.6 

± 

0.9 

± 

0.8 

Data 



107 





54 





6 



p-value 



0.36 





0.12 





0.24 






SRVLQ3 



SRVLQ4 


tmjz 

2.07 

± 

0.10 

± 

0.89 

3.14 

± 

0.13 

± 

1.35 

tin 

0.40 

± 

0.04 

± 

0.07 

0.57 

± 

0.05 

± 

0.07 

Dibosons 

2.36 

± 

0.29 

± 

0.61 

2.03 

± 

0.25 

± 

0.49 

Fake/Non-prompt 

3.09 

± 

1.29 

± 

1.80 

4.24 

± 

1.59 

± 

2.47 

Q mis-Id 

1.72 

± 

0.22 

± 

0.63 

1.45 

± 

0.17 

± 

0.52 

Other bkg. 

0.22 

± 

0.08 

± 

0.03 

0.41 

± 

0.10 

± 

0.06 

Total bkg. 

9.87 

± 

1.35 

± 

2.10 

11.9 

± 

1.6 

± 

2.8 

Data 



7 





10 



p-value 



0.83 





0.71 




Table 8 : Observed and expected numbers of events with statistical (first) and systematic 
(second) uncertainties for five of the signal regions defined for VLQ, chiral 5'-quark and 
four-top-quark production searches. The p-values for agreement between the observed yield 
and the expected background in each signal region are reported. 


CLfo are defined as the fractions of signal plus background and background-only pseudo¬ 
experiments with Lr larger than the observed value. Signal cross sections for which 
CLg = CLs+b/CLfc < 0.05 are deemed excluded at the 95% CL. Expected limits assum¬ 
ing the absence of signal are also computed; these are the basis for assessing the intrinsic 
sensitivity of the analysis. 

In the signal regions defined for searching for positively charged top quark pair produc¬ 
tion, the observed yields agree well with the expectation from background. The resulting 
limits on the cross section for this process are shown in table 10 in both the contact inter¬ 
action and Higgs-like FCNC models. For the special case of the 125 GeV Higgs boson, the 
limit on the cross section leads to a limit of BR(t —)• uH) < 0.01. The results can also be 
expressed as limits on the parameters defined in equations 1.2 and 1.3: for each chirality, 
the upper limit on C as a function of A is shown in figure 4; the same figure also shows the 
limits on HutH and KctH in the Higgs-like FCNC model. 

In contrast to the same-sign top signal regions, some of the signal regions defined for 
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SRVLQ5/SR4t2 

SRVLQ6/SR4t3 

SRVLQ7/SR4t4 

ttw/z 

1.87 

± 

0.09 

± 

0.80 

2.46 

± 

0.11 

± 

1.06 

0.57 

± 

0.05 

± 

0.25 

tlH 

0.31 

± 

0.04 

± 

0.05 

0.44 

± 

0.04 

± 

0.06 

0.08 

± 

0.02 

± 

0.02 

Dibosons 

0.33 

± 

0.14 

± 

0.10 

0.04 

± 

0.12 

± 

0.03 

0.00 

± 

0.12 

± 

0.00 

Fake/Non-prompt 

1.03 

± 

0.97 

± 

0.60 

0.00 

± 

1.02 

± 

0.28 

0.04 

± 

0.83 

± 

0.24 

Q mis-Id 

1.17 

± 

0.16 

± 

0.38 

1.09 

± 

0.14 

± 

0.34 

0.30 

± 

0.09 

± 

0.10 

Other bkg. 

0.16 

± 

0.08 

± 

0.02 

0.23 

± 

0.08 

± 

0.05 

0.14 

± 

0.08 

± 

0.08 

Total bkg. 

4.9 

± 

1.0 

± 

1.0 

4.3 

± 

1.1 

± 

1.1 

1.1 

± 

0.9 

± 

0.4 

Data 



6 





12 





6 



p-value 



0.46 




0.029 



0.036 



Table 9: Observed and expected numbers of events with statistical (first) and systematic 
(second) uncertainties for three of the signal regions defined for VLQ, chiral 6 '-quark and 
four-top-quark production searches. The p-values for agreement between the observed yield 
and the expected background in each signal region are reported. 


Model 

cr{pp —>■ tt) [fb] 

Coupling const. 


Exp. 

Obs. 

Observed 

Contact interaction model 

|C|/A2 [TeV-2] 

Left-left 

64 

62 

0.053 

Left-right 

53 

51 

0.137 

Right-right 

40 

38 

0.042 

Higgs-like FCNC model 

KutH or KctH 

uu —>■ tt [mu = 125 GeV) 

37 

35 

0.16 

uu —>■ tt [mu = 250 GeV) 

21 

20 

0.17 

uu —>■ tt {niH = 500 GeV) 

12 

11 

0.20 

uu tt {rriH = 750 GeV) 

9.3 

8.4 

0.24 

cc tt {niH = 250 GeV) 

71 

69 

0.81 

cc tt [mn = 500 GeV) 

37 

35 

1.02 

cc ^ tt {niH = 750 GeV) 

28 

27 

1.29 


Table 10: Observed and expected 95% CL upper limits on the cross section for same-sign 
top-quark production, and on the coupling constants. 


VLQ, 5'-quark, and four-top-quark production exhibit an excess over expected background. 
The excess is largest in the subset of the signal regions used for the four-top-quark search, 
where at least two 6 -tagged jets are required. While it is still of interest to limit the set 
of models consistent with the data as described above, it is also important in this case 
to assess the consistency of the data with the background-only hypothesis. This is done 
by computing p = 1 — CLf,. The resulting p-values depend on the signal model and the 
signal regions considered, as shown in figures 5a and 5b. For signals where all eight signal 
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Figure 3: Expected background yields and observed data events in various signal regions. 
Uncertainties include both the statistical and systematic errors. The difference between 
data and expectations is quantihed by the means of the significance, computed from the 
p-values in tables 8 and 9. 


regions are considered (as is the case for VLQ and V models), the significance is above 
one standard deviation but less than two. For signals for which only SR4tO-SR4t4 are 
considered (as is the case for four-top-quark production models) the significance reaches 
2.5 standard deviations. Several checks (detailed in section 8) of the background estimates 
were performed. Some features of the events in the signal regions that exhibit the most 
significant excesses (SR4t3 and SR4t4) are presented in section 9. 

The excess is not significant enough to support a claim of BSM physics. Therefore 
95% CL limits (upper limits on cross sections, or lower limits on masses) relevant for each 
model are calculated. The observed excess causes these limits to be less restrictive than 
expected for the background-only hypothesis. The data place 95% CL upper limits on 
the 5'-quark pair production cross section that vary with the mass of the 6'-quark. Limits 
obtained assuming a 100% branching ratio to Wt are presented in hgure 6a (expected mass 
limit at 0.79 TeV, observed at 0.73 TeV), and limits where decays to u- or c-quarks are also 
considered are shown in figures 6b and 6c. 

Limits on the VLQ pair-production cross section, assuming the branching fractions to 
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Figure 4: Expected and observed limits on the coupling constant \C\ in the contact 
interaction model for same-sign top quark pair production as a function of the BSM physics 
energy scale A for the (a) left-left, (b) left-right and (c) right-right terms. Plot (d) shows 
the observed limits on the Higgs-like exchange with FCNC coupling in the plane {k-uih■,K- ctH) 
for three different hypotheses on the mass of the heavy Higgs-like particle. 


VF, Z, and H modes prescribed by the singlet model, are shown in figure 7. Comparison 
with the calculated cross-section results in lower limits on the H-quark mass of 0.62 TeV 
and on the T quark mass of 0.59 TeV at 95% CL. The expected limits in the absence of 
a signal contribution are 0.69 TeV for the H-quark mass and 0.66 TeV for the T-quark 
mass. If the three branching fractions are allowed to vary independently (subject to the 
constraint that they sum to one), the data can be interpreted as excluding at 95% CL some 
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(a) (b) 

Figure 5: (a) Probability for the data in the four-top-quark signal regions (SR4tO-SR4t4) 
to be consistent with a zero cross section for anomalous four-top-quark production under 
two model scenarios, as a function of the characteristic mass scale of the models, (b) 
Probability for the data in the VLQ/ 6 '-quark signal regions (SRVLQ0-SRVLQ7) to be 
consistent with a zero cross section for various heavy quarks, as a function of the quark 
mass. 


of the possible sets of branching ratios for a given B- or T-quark mass. These exclusions 
are shown in figures 8 and 9. 

Limits on T 5/3 production are set for pair production only, and for the sum of pair 
and single production for two different values of the coupling A of the T 5/3 to Wt (A = 0.5 
and 1.0) [89]. This coupling is related to the mixing parameter g* used by the model in 
refs. [90, 91]: A = niT^^^gg* j mwV^- The pair-production limits are shown in figure 10a, 
and correspond to a mass limit of 0.74 TeV (0.81 TeV expected). The limits on pair plus 
single production with A = 0.5 are shown in figure 10b, where the observed mass limit is 
0.75 TeV and the expected limit is 0.81 TeV. Finally, limits on pair plus single production 
with A = 1.0 are shown in figure 10c, where again the observed mass limit is 0.75 TeV and 
the expected limit is 0.81 TeV. 

The upper limit on the cross section for four-top-quark production is 70 fb assuming 
SM kinematics, and 61 fb for production with a BSM-physics contact interaction (expected 
limits are respectively 27 fb and 22 fb). The cross-section limit for the contact interaction 
case is lower than for the SM since the contact interaction tends to result in final-state 
objects with larger pxj which increases the selection efficiency. The limits are also inter¬ 
preted in the context of specific BSM physics models. For the contact interaction model, 
the upper limit on |C' 4 t|/A^ is 15.1 TeV“^, as illustrated in figure 11a. The lower limit on 
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Figure 6: Limits on ft'-q^ark pair production: (a) Expected and observed upper limits 
on the h' pair-production cross section times the branching ratio for b'U as a 

function of the 6'-qaark mass (the vertical dashed lines indicate the expected and observed 
limits on the 5'-qaark mass, and the shaded band around the theory cross section indicates 
the total uncertainty on the calculation); (b) expected and (c) observed exclusion limits on 
the 5'-qaark mass as a function of the assumed branching ratios into c- and t-quarks. 


the sgluon mass is 0.83 TeV, assuming that the sgluons are pair-produced and always decay 
to ti (for an expected limit of 0.94 TeV), as shown in figure 11b. The observed limits on the 
cross section times branching ratio for the 2UED/RPP signal are shown in figures 11c, lid 
and 12. These imply the following limits on tukk: in the symmetric case (i ?4 = R 5 ), the ob¬ 
served limit coming from tier (1,1) is 0.96 TeV (where the expected limit is 1.05 TeV). The 
observed limit coming from tiers (2,0) -|- (0, 2) alone (BR(Ai^d) |_^ = g) jg o_ 5 q TgV 
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(a) 


(b) 


Figure 7 : Expected and observed limits on the pair production cross section as a function 
of mass for (a) vector-like B and (b) vector-like T quarks. The vertical dashed lines 
indicate the expected and observed limits on the vector-like quark mass. These limits 
assume branching ratios given by the model where the B and T quarks exist as singlets [25]. 
The shaded band around the theory cross section indicates the total uncertainty on the 
calculation. 



Figure 8: Expected (a) and observed (b) vector-like B quark mass hypotheses excluded 
at 95% CL as a function of the assumed branching ratios. 


(where the expected limit is 0.55 TeV). In the highly asymmetric case (i ?4 > R 5 ), tier 
(0, 2) does not contribute any longer and the observed limit on ut-kk from tier (2, 0) alone 
is 0.45 TeV (where the expected limit is 0.51 TeV). Eigure 12 shows the limits in the mKK~C 
plane, with the constraints from cosmological considerations superimposed. 
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Figure 9: Expected (a) and observed (b) vector-like T quark mass hypotheses excluded 
at 95% CL as a function of the assumed branching ratios. 
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Figure 10: Expected and observed cross section limits as a function of mass on (a) 
pair production, (b) T 5/3 pair plus single production for coupling A = 0.5, and (c) T 5/3 
pair plus single production for A = 1.0. The vertical dashed lines indicate the expected 
and observed limits on the T 5/3 mass, and the shaded band around the theory cross section 
indicates the total uncertainty on the calculation. 
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Figure 11: Limits obtained from the search for four-top-quark production, (a) Expected 
and observed limits on the coupling constant | C^t \ in the contact interaction model for four- 
top production as a function of the BSM physics energy scale A. The region in the top left 
corner, corresponding to |C' 4 t|/A^ > 15.1 TeV“^, is excluded at 95% CL. (b) Expected and 
observed limits for the sgluon pair-production cross section times branching ratio to four 
top quarks as a function of the sgluon mass, (c) Expected and observed limits on the four- 
top-quark production rate for the 2UED/RPP model in the symmetric case. The theory 
line corresponds to the production of four-top-quark events by tier (1,1) with a branching 
ratio of ttti of 100%. (d) Expected and observed limits on the four-top-quark 

production rate for the 2UED/RPP model in the symmetric case. The notation (2) -|- (2) 
is a shorthand for -I- Aib’ 2 )j 4 {o, 2 )^ theory line corresponds to the four-top- 

quark production by tiers (2,0) -|- (0, 2) alone (BR(Ai^T) = q) q;’]^g vertical dashed 

lines in (b), (c), and (d) indicate the expected and observed limits on the sgluon mass or on 
rriKK, and the shaded band around the theory cross section indicates the total uncertainty 
on the calculation. 
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Figure 12: Expected and observed limits in the (tukk = 1 /-R4,^ = Ra/R^) plane for the 
2UED/RPP model. Cosmological constraints are also shown. They are due both to direct 
searches for dark matter and from constraints on its relic density (see ref. [43] for more 
details). The mioc parameter from ref. [43] affects the dark-matter relic density predicted 
by the 2UED/RPP model. 
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8 Checks of the background estimate 


Several checks were performed to assess the validity of the background estimate. The most 
important tests are summarized here. For the simulation-based background estimates: 
variations in the cross section, in the generators, and their settings that span the range 
consistent with theoretical expectations or direct measurements were applied. Variations 
in the expected yield are part of the systematic uncertainty. The data-driven estimates 
(for the charge misidentification and fake/non-prompt lepton backgrounds) were checked in 
several ways. The particular leptons observed in SRVLQ6 and SRVLQ7 were scrutinized, 
and their quality was found to be consistent with a sample dominated by real leptons. Sim¬ 
ilarly, the multivariate discriminant for 5-tagging is well above the required threshold for 
tagged jets found in the sample. In addition, the expected contribution from charge misid¬ 
entification and from fake/non-prompt leptons was assessed using samples of simulated 
events. It is found that the yields are consistent within uncertainties with the expectations 
from the data-driven estimates. To further investigate whether the matrix method accur¬ 
ately predicts the number of fake/non-prompt leptons in tt events, the entire procedure 
was repeated with samples of simulated events, with r and / measured using simulated 
single-lepton and multijet samples respectively. The predicted number of fake/non-prompt 
leptons in simulated tt samples is consistent with the actual number present in the MC 
samples. 

9 Features of events in signal regions with most significant excesses 

Information about the lepton charges and flavours, as well as some key kinematic inform¬ 
ation, for the observed events in signal regions SR4t3/SRVLQ6 and SR4t4/SRVLQ7 are 
presented in tables 11-12. One unexpected feature is the dominance of one electric charge 
over another: in SR4t3/SRVLQ6 there are 10 negatively charged and 16 positively charged 
leptons, and in SR4t4/SRVLQ7 there are only 2 negatively charged leptons and 11 posit¬ 
ively charged leptons. Similar effects are observed in two other signal regions, with SRVLQ3 
being dominated by negatively charged leptons and SR4t2/SRVLQ5 being dominated by 
positively charged leptons. These charge asymmetries are interpreted as statistical fluctu¬ 
ations. This interpretation is bolstered by the fact that i) there is no known mechanism 
for the selection to favour one electric charge over the other; ii) the asymmetry is present 
only in some of the signal regions (and both negative and positive charges dominate in the 
various regions); and iii) that the asymmetries do not persist in the loose lepton samples 
selected with the same kinematic criteria as are applied in the signal regions. 
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Type 


Ht [GeV] 

[GeV] 

e~ e~ 

3 

807 

171 

e+e+ 

5 

862 

268 

e+e+ 

5 

868 

113 

li~e~ 

6 

1346 

353 


5 

810 

106 

e~p~ 

3 

707 

184 

e~p~ 

2 

706 

174 


8 

882 

150 

^+e+ 

4 

860 

112 


5 

888 

111 

pL~e'^e^ 

5 

773 

197 

pL~e'^e^ 

9 

968 

355 


Table 11: List of data events in the SR4t3/SRVLQ6 category (by definition, all events 
have exactly two 6-jets, the is above 100 GeV and the above 700 GeV). 


Type 

Nj 

Nb 

Ht [GeV] 

[GeV] 

e+e+ 

4 

3 

709 

298 

e+e+ 

6 

3 

800 

137 

e+^+ 

5 

3 

744 

216 

e+;U+ 

4 

3 

888 

155 

;U+e+ 

3 

3 

1439 

239 


4 

4 

1072 

176 


Table 12: List of data events in the SR4t4/SRVLQ7 category (by definition, the is 
above 40 GeV and the TLt above 700 GeV). 

10 Conclusion 

A search for BSM physics has been performed using pp collisions at = 8 TeV corres¬ 
ponding to an integrated luminosity of 20.3 fb“^ recorded by the ATLAS detector at the 
LHC, where events with at least two leptons, including a pair of the same electric charge, 
at least one 6-tagged jet, sizeable missing transverse momentum, and large 77 t were con¬ 
sidered. Several BSM physics effects could enhance the yield of such events over the small 
SM expectation. The search was performed in the context of several BSM physics models, 
with signal regions defined for different models. The regions of parameter space excluded 
by the data are quantified by setting 95% CL limits. The observed yield in the signal 
region for positively charged top quark pair production is consistent with the expected 
background, resulting in limits of 8.4-62 fb on the cross section of this process (depending 
on the model considered) and a limit BR(t —)• uH) < 1%. In the set of signal regions 
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defined for vector-like quark, four-top-quark, and chiral 6 ^-quark searches there is an excess 
of observed events over the SM prediction, particularly in the subset of those signal regions 
that require at least two 6 -tagged jets (and thus are relevant to the search for four-top- 
quark production). The significance of the excess varies with the signal being considered, 
reaching 2.5 standard deviations for hypotheses involving heavy resonances decaying to 
four top quarks. Nonetheless the data can still constrain some of the BSM physics models 
considered; 95% CL limits are set as follows: the mass of the chiral 6 Lquark is constrained 
as a function of the branching ratio to Wt, the masses of vector-like B and T quarks are 
constrained to mB > 0.62 TeV, rriT > 0.59 TeV (assuming branching fractions to the W, 
Z, and H decay modes arising from a singlet model), the mass of the T 5/3 quark is greater 
than 0.75 TeV, the SM four-top production cross section is less than 70 fb, the sgluon mass 
is greater than 0.83 TeV and the Kaluza-Klein mass (in the context of models with two 
universal extra dimensions) is greater than 0.96 TeV. 
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